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Abstract 

By solving the gap equation in the quenched, ladder approximation for an 
Abelian gauge model with Yukawa interaction in the presence of a constant 
magnetic field, we show that the Yukawa interactions enhance the dynamical 
generation of fermion mass. The theory is then studied at finite temperature, 
where we prove that the critical magnetic field, required for the mass genera- 
tion to be important at temperatures comparable to the electroweak critical 
temperature, can be substantially decreased due to the Yukawa coupling. 

I. INTRODUCTION 

In recent years the phenomenon of dynamical symmetry breaking in the presence of an 
external magnetic field has attracted a great deal of attention |[|- [[|. The essence of this 
effect lies in the dimensional reduction in the dynamics of fermion pairing in the presence 
of a magnetic field [jj. Due to such a dimensional reduction, the magnetic field catalyses 
the dynamical generation of a fermion condensate and a fermion mass, even in the weakest 
attractive interaction between fermions. 

The field-induced dynamical generation of mass (FDGM) has been found in many models 
of field theories, in 2+1 and in 3+1 dimensions. The universality of this phenomenon makes 
it of interest for a wide set of applications, including problems in different areas such as 



cosmology ||, [[|, astrophysics JT0|, and planar condensed matter ]TT[ |T2 . 



In their original papers [|1], Gusynin , Miransky and Shovkovy suggested that the struc- 
ture of the electroweak phase transition could be affected by the FDGM. Considering FDGM 
in QED 4 at finite temperature, Lee, Leung and Ng JJ, and independently, Gusynin and 
Shovkovy ||, calculated the critical temperature at which the fermion mass (and hence the 
fermion condensate) evaporates. From their results one might conclude that the dynamical 
generation of mass induced by a magnetic field should play no role in the electroweak phase 
transition, because for it to exist at the high temperatures typical of the electroweak scale, it 
would require a primordial magnetic field too large (~ 10 42 G) to be realistically attainable. 
However, as we have argued in a previous paper ||, it is reasonable to expect an essential 
modification in the order of the critical field when the FDGM takes place in the context 
of the electroweak model, since a richer set of interactions enters in scene there. As we 
will show below, this is indeed the case even in a model much simpler than the electroweak 
theory. By studying a toy model similar to QED 4 , but including a Yukawa term, we prove 
that the Yukawa interactions can substantially enhance the FDGM phenomenon. 

Any cosmological application of the FDGM has to assume that primordial magnetic fields 
could have been present during the early times of universe evolution. Nowadays astronomical 
observations seems to support this assumption. As it is known, primordial magnetic fields 
may be needed to explain the large-scale galactic magnetic fields ~ 1CT 6 G observed in 
our own, as well as in other galaxies. The observed galactic fields have been a source 



of motivation for many works on primordial- field generating mechanisms ||13|| . Typically, 
fields as large as 10 24 G are predicted by these mechanisms during the electroweak phase 



transition. Moreover, Ambjorn and Olesen jrjj have claimed that seed primordial fields 
even larger, ~ 10 33 G, would be necessary at the electroweak scale to explain the observed 
galactic fields. 

From the above discussion it can be understood that the significance of the FDGM at the 
electroweak scale needs yet to be elucidated. The present paper is a step in that direction. 
We study the FDGM in an Abelian gauge model of massless fermions with a Yukawa term. 



This model, although simple, retains some of the attributes of the electroweak theory. We 
explicitly show that Yukawa interactions enhance the generation of mass in the presence 
of a magnetic field, decreasing the critical field needed for the FDGM to be important at 
temperatures comparable to the electroweak critical temperature. For a Yukawa coupling 
of order of the top-quark coupling, the critical field strength is decreased in 10 orders of 
magnitude as compared to the corresponding field strength in QED. 

The plan of the paper is the following. In Sec. II, we derive the Schwinger-Dyson 
equation for the fermion self-energy in the quenched, ladder approximation, solving the 
corresponding gap equation in the infrared region. Thermal effects are considered in Sec. 
Ill, where we calculate the critical temperature at which the field-induced fermion mass 
disappears and estimate the order of the critical field required for the FDGM to be significant 
at the electroweak scale. We discuss the implications of our results and state our conclusions 
in Sec. IV. The solution of the SD equation for the coefficient Z\\ of the mass operator is 
derived in the Appendix. 



II. ABELIAN GAUGE MODEL WITH YUKAWA INTERACTION 

Let us consider an Abelian gauge model with a Yukawa interaction described by La- 
grangian density 

L = --F^F^ + - #r"V>A* - \^A^f + ~d M <W - ^ 4 - V2\ y m (1) 

Note that L has a U(l) gauge symmetry and a fermion number global symmetry, but 
it does not have a continuous chiral symmetry. We consider the present model in the 
presence of a constant and uniform external magnetic field H . Our aim is to investigate 
the gap equation of the theory. We need to go beyond a perturbative calculation to find 
a non zero solution of the gap equation. In the small coupling regime of the theory a 
consistent, non perturbative calculation can be carried out by using the quenched, ladder 
approximation. This approximation has produced gauge invariant \TE\ non trivial solutions 
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of the gap equation in other theories with external fields so it is natural to expect that 
it will yield similar results in the present case. 

The gap equation of the theory described by the Lagrangian density (|lj) can be obtained 
from the Schwinger-Dyson equation for the fermion self-energy 

G 1 (x,y) = G' 1 {x,y) +ig J d A ud A wyG(x,u)D(x - w)T^ A (u,y,w) 

+ iV2\ y J d 4 ud A wG(x,u)S(x - w)T^(u,y,w) (2) 

Here G refers to fermion propagators, and D and S to gauge and scalar boson propagators 
respectively. and are three-fields vertex functions. The bar indicates full Green 

functions. A compact notation where tensorial and spinorial indexes have been suppressed 
is understood. 

After taking the quenched, ladder approximation (where the fermion propagator in the 
presence of an external magnetic field is taken full, while the vertices, as well as the gauge 
and scalar boson propagators, are taken bare) of Eq. (Q), one obtains the following equation 
for the fermion mass operator 



M(x, y) = G (x, y) — G (x, y) = —ig J d ud wyG(x, u)^D(x — w) 

-i (V2\ y ) 2 J d 4 ud A wG(x, u)S(x - w) (3) 

To solve Eq. (|3|) we need to transform it to momentum space. However, it is known 
that the Fourier transform of the fermion Green function in the presence of a magnetic field 
does not yield a diagonal-in-p function. The reason is that in the problem with external 
field the fermion asymptotic states are no longer plane waves. A suitable solution to this 



technical problem was found many years ago by Schwinger [RJ, who introduced the proper- 
time representation of the fermion Green function in the presence of a constant external field. 
Here, however, we prefer to adopt another approach due to Ritus |L7[]. Ritus' method is based 
on the use of a representation (known as the E p representation) spanned by the solutions ip p 
of the eigenvalue equation (7n) 2 ^ p = p 2 ip p . The operator (7n) 2 = (7 M (z<9 M — gA /t )) 2 commutes 



with the mass operator. Hence, the ipp can be used to generate a set of complete and 
orthonormal eigenfunction-matrices of the mass operator. In the chiral representation of the 
Dirac matrices, where 75 and ^3 = 27172 are both diagonal, the ifi p take the form 

ip P = E pax (x)u ax (4) 

The bispinors uo ax are eigenvectors of 75 and with eigenvalues x = ±1 and a = ±1, 
respectively. 

In the case of a purely magnetic field background (crossed field case) directed along the 
z-direction, the E pcrx functions are 

E pa (x) = N{n)e l(pox ° +P2x2+P3x3) D n {p) (5) 



where D n (p) are the parabolic cylinder functions [|T1| with argument p = y2 \gH\(x\ — jfj] 



and positive integer index 

gHa 1 

n = n(k,a) = k + - — — - -, n = 0,1,2,... (6) 

N(n) = (An \gH\)^ /y/nl is a normalization factor. Here p represents the set (po,P2,P3, k), 

which determines the eigenvalue p 2 = —p\ + p\ + 2 \gH\ k in ('-yTVf'ipp = p 2 ip p . Note that in 
this case E pax does not depend on x- 

The E p representation is obtained forming the eigenfunction-matrices 

E p (x) = Y J E p(T (x)A(a) 1 (7) 

(7 

where 

A(cr) = diag(8 al , S al , ° = ±1, (8) 

It is easy to check that the E p functions are orthonormal 

' d A xE pl (x)E p (x) = (2n)^ 4 \p - p') = (2n)% k ,5(p - p' )5(p 2 - p' 2 )5(p 3 - p' 3 ) (9) 
as well as complete 



X; / dp dp 2 dp 3 E p (x)E p (y) = (2n) 4 5^(x - y) (10) 
k J 

Here we have used E p (x) = ■y°E^'j . 

They also satisfy two important relations 

7 ■ UE p {x) = E p (x)i ■ p (11) 

J d 4 x'M(x, x')E p {x) = E p (x)Ea(p) (12) 

where Y^a{p) is the fermion mass operator in momentum coordinates. 

In the E p representation the SD equation (||D for the mass operator becomes 

(2vr) 4 (W<5(po - Po)$(P2 - p' 2 )S(p3 - p' 3 )Ea(p) 

dp\dp\dp" 3 r - 



V/ d^x'Y. J °Zw 3 {e p (^eax 



(2tt) 4 PV 7"P"-Ea(P") 



x^^f^M^^ - *')} - z2Aj j AdV^ | dfo ^ f " {E p (x)E lf (x) 



x _„ ~ j— EfWEjWSix - x')} (13) 



with 



D,u(x -x') = -[ ( d X e ~T' X ' ] L» - (1 - O^A (14) 
the bare photon propagator, and 

f d 4 q e-^-O-z') 
S ^-x)=j (27r) 4 g2 _, e (15) 

the bare scalar propagator. 

Using the properties of the parabolic cylinder functions and Eqs. (|TT|)-(|T2|), the integrals 
in x and x' in Eq. (|T3|), as well as the integrals in p 0j p 2 ,&nd p 3j can be done yielding 



~ f d^n p isgn(gH)(n-n"+n" -n')tp i 



Ml) a 7 ^a" 1 

? I 7-P"-Ea(P" 



x ( g, v - (1 - ] A 7 ^A" - _ ^ A VA'} 



(i 4 (7 e isgn(gH)(n-n n +n" -n')tp ^ 

e" 9 -L Jnn"(q± ) 



xAA" = A" A'} (16) 

7"P"-Ea(P") 

where the following notation has been used 

min(n p ,n r ) | | 

Jn P nM±)= E — ( ^9n(9H)q ± r^- 2m (17) 

^ m\(n p — m)\(n r — my. 



p" = (po - g , 0, -sgn{gH)p \gH\ k",p 3 - q 3 ) (18) 

E= E^ (iQ) 

{°"} CTCr"(T / (T" 

n = n(Jfe,<r), n' = n(k',a'), n"=n{k",a"), n" = n(Jfe" ,?"), 
and the dimensionless variables 



? M = V 2 h ' A* = 0,1,2,3 (20) 



and polar coordinates for the transverse components of g M 



9i = V§i + V 9 = arctan(g 2 /gi) (21) 

have been introduced. 

The first term of the right hand side (RHS) in Eq. ([16]) coincides with the result previ- 
ously found in references [3] for the pure QED case. Because of the presence of the Yukawa 
interaction in the present model, we have obtained an additional contribution (the second 
integral in the RHS of Eq. flTBD) whose consequences for the gap equation will not be trivial, 
as we show below. 



7 



To solve equation (0), we need the mass operator structure. In the presence of the 
external magnetic field the mass operator structure is quite rich fl5(| . However, as argued 
in a previous paper [[yj, within the present approximation a simpler structure can be used. 
Thus, we consider 



(22) 



Note the separation between transverse and perpendicular variables. 

After substituting the above structure for J2a{p)i Eq. (|T6| ) can be written as 



Zf(-p + Z ± >y-p ± +m( 



ig 2 2\gH\^[ 



e isgn(gH){n-n ,, +n' , -n')ip 



frfaJ (2tt) 4 Vn\n n \n n \n'\ 



c o r . 



{( aa\ - (l + Z u h-p" - (l + ZA^-p" _ 



1 + Z ) 2 p" 2 + (1 + Z ± ) 2 p 2 + m 2 (p") 



-.2A;(2|^|)5:E / 



d A q e is 9n(9H)(n-n"+n"-n')<fi ^ ^ 
^ / /^I / , ,,,-,,, =77 e 9x ^nn"(§±)^n"n'(9J-)^ 

jftriJ (27TJ 4 i/n!n !n !n'! <T 



f m(p" ) - (1 + Z ) 7 • p " - (1 + Z ± ) 7 • p " „ ; 
x AA" — : — -1— ; — i — : — — ^— A "A 



(23) 



(1 + Z v ) 2 p" 2 + (1 + Z L fp\ + m 2 (p") 
It can be further simplified taking into account that large q± contributions to Eq. (|23|) are 
suppressed by the factor e - ^, and therefore we can approximate the J nn '(q±) functions by 
their infrared behavior 



J nn '(q±) ~ n\5„ 



(24) 



This approximation allows us to eliminate the tp dependence in the integrand to obtain 



<W Zff ■ P., + Z17 • p ± + m{ 



zg 2 2\gH\Y,Y. S nn»k» n >J 



d A q e 9 -l 



k" {a} 



(2vr) 4 g 2 V(1 + ^|,) 2 P" 2 + (1 + ^±) 2 P 2 +m 2 (j>"' 



x{m(p n 

-a+^i,) 
-o-+z±: 



AyA'A'fW - ^^ A( 7 • g)A"A"( 7 • q)A' 



A 7 ^A" ( 7 .p 1 »)A"7 /l A 



, (i-0 



<7 2 



A( 7 -g)A" ( 7 -p || ")A"( 7 -g)A / 



A^A" ( 7 -p ± ")A" 7m A 



, (1-0 



g 2 



A( 7 -g)A" ( 7 -p ± ")A"( 7 -g)A' 



} 



-i2Aj(2|s#|)££(W(W / 



fc» {a} 



d 4 g e 9 -l 
(2vr) 4 g 2 



x 



m(p")AA"A"A / - (1 + Z„)AA" ( 7 • A" A' - (1 + Z ± )AA" ( 7 ■ p ± ") A"A /X 



(25) 



(1 + Z^p" 2 + (1 + Z ± ) 2 p 2 ± + m 2 (p") 

To perform the summation in the spin indices we can make use of the following relations 
satisfied by the A matrices (H) 



A(a)A(a') = 5^A{a) 



(26) 



A(l) + A(-1) =/ 



At? = 7f A, 



A 7 ^ = ^(1 -A), 



with 



7 f = ( 7 °,0,0, 7 3 ) 



tyitft ^ = (0,^,^,0) 



(27) 
(28) 



to find that 



E 

M 



]T 5 nn »<^>, n ,AA"A"A' = <W<W (29) 



(W^A^A" ( 7 • p ± ») A" 7m A' = 2 ( 7 • <W<W (30) 



^<W<W A 7 M A" ( 7 • p ") A" 7(U A' = 2 ( 7 • p ») 8 W (5 kr _ sgn{gI1) A{l) + 5 kr+sgn{gH) A(-l 
M 

(31) 
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J2 <W«WA ( 7 • q) A" (7 • p±" ) A" (7 • g) A' ~ (7 • p ± ") <W<W 



(32) 



A( 7 -g)A" (7-P„")A" ( 7 -g)A' (7 • g„) (7 • P„») (7 ■ ^) o 

2^ d nn"^"„' ^ — <Jkk><Jkk" ^ [OO) 



X; <W<%„ n ,AA" (7 • p x ") A" A' = <W<W (7 • 
M 



(34) 



E^-W^" (7 • v) A " A ' = 4fc,4fc " ( 7 -Pii 



(35) 



In the above equations we dropped terms proportional to q ± taking into account the 
small ^approximation here considered. The remaining a sums were obtained in ref. |§. 
Note that the 8^1 appearing in all terms cancels out with the one of the left hand side 
(LHS) of Eq. fl25l) . The SD equation is then expressed as 

d 4 q ej / 1 \ 

(2tt) 4 q 2 + Z..) 2 f' 2 + (1 + Z L ) 2 p 2 + m 2 (p") J 



Z ][ j-p ii +Z±-fp ± +m(p) = i2\gH\J2 J 



X{g 2 {m(f) (1 - f) 5 fefc " - 2 + 5 k ,k"-sgn(gH)A(l) + Sk,k»+sgn(gH) A(-l 

(7-9,,) (t-P,,") (7-?j|) 



+(1 + z„ 



(i-O^fc- 



— 2 (7 •£>,,") (£jfc,fc"-sgn(g-ff)A(l) + <5jfc,fc'>+ sgn (gff)A( 



+ (1 + Z ± ) [(1 - <W (7 ■ P±n - 24 fe " (7 • P ± ")] } 

-2Aj{m(?f )<W - (1 + Z„)<W (7 • £„") - (1 + Z ± )5 W , (7 ■ P ± )}} 

The A's terms, appearing in the RHS of Eq. (^), but not in the LHS, remind us that 
the mass operator structure used here, Eq. (|22|), is not the most general one in the presence 



of an external magnetic field fl5|| . However, to study the dynamical generation of mass we 
can restrict our calculations to certain momentum region on which the simplified structure 
([TT|) gives rise to self-consistent results. Let us recall that in the case of QED it has been 
shown [|TJ , , p| that the dynamical generation of mass in the presence of a magnetic field 
is governed by the fermion infrared dynamics. It is natural to expect that a similar situation 
takes place in the model we are considering here. Thereby, we shall restrict the above SD 
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equation to the lower Landau level (LLL) approximation, on which p = k = 0, which is 
better justified if we simultaneously assume that the external momentum lies in the infrared 
region p 2 « \eH\ . If, in addition, we use the Feynman gauge, and realize that the main 
contribution to the sum in k" comes from the k" = k = terms, we can rewrite the SD 
equation as 

Z {{1 ■ p B + m(p„) - ^2 \ 9 H\ / ( (1 + Z )^ + m^p" ) ) 



x 



{-2(g 2 + Aj)m(jf ) + 2A}(1 + ^)<W ( 7 -p,,")} (37) 



Performing a Wick rotation to Euclidean space, Eq.fl37|) leads to the following two equa- 
tions, one for each independent structure, 

r d 4 d p-Q± (1 + Z.M ■ (p a — q ) 

°.t*, - -«/(^V i(i +w ,\ )a : m 4-^i (38) 



m(p, ) = 4 |</if| (<? 2 + X 2 ) / — — T (39) 



(2tt)V ( i + Z)2(p )2 + m 2 (p 



As shown in the Appendix, the first equation has solution Z. = 0. The second is the gap 
equation. Note that it is very similar to the gap equation found by Lee, Leung and Ng || 
(after using Z. = 0), except that the coupling factor here is (g 2 + X 2 ) instead of just g 2 . The 
main contributions to Eq. (|39|) come from the infrared region q\ « \gH\ , q 2 << \gH\ . 
A solution to Eq. ([59]) can be explicitly found in the infrared approximation p ^ 0, taking 
into account that the mass parameter is dominated by the small momenta contributions || . 
Then, we can approximate the mass by its infrared value m(p^ — q ) ~ m(p^) ~ m(0), and 
use that Z.. = to arrive at the infrared gap equation 

\qH\(q 2 + X 2 ) 7 7 e -%- e~*± 

1 ~ 19 l[g + y> / d 2 a i / d 2 a — (40) 

An 2 J q± J q wql + ^2\gH\q 2 + m 2 ' 1 ' 

which can be integrated to obtain the non trivial solution 
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m ~ 



\J\gH\exp 




7i 



(41) 



,2 



The consistency of the approximation requires m « J\gH\, which is satisfied if + 

A 2 

<< 1, so the dynamical mass appears in the weak coupling region of the theory. Note 
that, because of the exponential function in Eq. fl4lf) , small changes in the exponent can 
yield substantial changes in the mass. For instance, for X y ~ 0.7, a value comparable to the 
top-quark Yukawa coupling, the dynamical mass is five orders of magnitude larger than 



We should underline that in this model the generation of a dynamical mass cannot be 
linked to the breaking of a continuous chiral symmetry, because there is no chiral symmetry 
to begin with. Since our goal is to use the present toy model to get insight of the FDGM 
phenomenon in the electroweak theory, where there is no chiral symmetry to break, the 
present model is a good candidate for our purpose. We must point out that, even though 
the appearance of the dynamical mass in the present case can be easily traced to the existence 
of a fermion-antifermion condensate there is no Goldstone field produced in this theory, 
because there is no continuous symmetry broken by the fermion condensate. We expect that 
this should not be the case in the electroweak model. If a fermion condensate is catalyzed 
by the magnetic field there, it might give rise to a field-dependent vev of the scalar field and 
hence to a Higgs-like spontaneous gauge symmetry breaking. 

Eq. ([|l]) clearly indicates that Yukawa interactions enhance the dynamical generation of 
the fermion mass in the presence of a magnetic field. To grasp the possible significance of 
this result in the context of the electroweak phase transition, finite temperature effects has 
to be considered. 



the mass found in QED [p]] — ||, which is given by m 




III. FINITE TEMPERATURE CALCULATIONS 



Finite temperature effects can be incorporated using the well known imaginary-time 
Matsubara formalism, in which the Euclidean time variable is compactified to a circle of 
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radius (3 = ^, (T is the absolute temperature), the fourth components of the momenta are 
consequently discretized according to 

g 4 = 2mi/P with n = 0, ±1, ±2, ... for bosons 

= (2n + 1)tt / /3 with n = 0, ±1, ±2, ... for fermions 

and, in the functional integrals, boson (fermion) fields are periodic ( ant iper iodic) in time 
with period (3. 

The calculation leading to the gap equation in the above section can be performed at 
finite temperature in a similar way, taking into account that now the integrals in the fourth 
components of the momenta must be substituted by sums due to the discrete character of 
these variables. Hence, the finite temperature gap equation can be expressed as 



m(u). 



n ' ,P) Mtt In'n^J ul + k 2 + m 2 (u n ,k)J (u n -u n ,y + (k- P y + x 1 } 

— oo 

with uo n = (2n+ l)ixT. Eq. ([4"2]) differs from the corresponding equation for QED (Eq.(10) of 
reference ||) only in the factor multiplying the integrals, which in the present case contains 

A 2 

the contribution ~r due to the Yukawa interaction. 

4-7T 

It is natural to expect that at some critical temperature the thermal effects evaporate the 
ifjip fermion condensate responsible for the nonzero dynamical mass. For the dynamical mass 
generation to be of any significance at the electroweak scale, it is needed that the critical 
temperature at which the mass becomes zero results of the order of the electroweak critical 
temperature ~ 10 2 Gev. Of course, we must remember that there is another parameter in the 
problem, the magnetic field. Different magnetic field strengths will yield different values of 
the critical temperature. Therefore, our goal is to determine the magnetic field required to 
have a critical temperature comparable to the electroweak one. Hence, we need to find the 
critical temperature as a function of the couplings and the magnetic field, that is, we need 
to solve first the gap equation ([4"2"D, and then take the limit m — > 0. Since our calculations 
are very similar to those of ref. |4j], we remit the interested lector to that paper, and here 

13 



we just give the final result, which in our case contains the modification due to the extra 
coupling X y . Thus, the critical temperature estimate is 



which results comparable to the dynamical mass at zero temperature. Critical temperature 
estimates of the order of the corresponding zero temperature dynamical masses have been 
also found in QED and NJL, in (2+1) and (3+1) dimensions Q § [0. 

We can now estimate the strength of the magnetic field required to have T c ~ 10 2 Gev . 



thus the critical field is H « 10 32 G. This result represents a decreasing of the field in 10 
orders of magnitude as compared to the value required in QED (where X y ~ 0, g = e) to 
obtain the same critical temperature of 10 2 Gev . 

In a note added in proof, Lee et. al. concluded that the dynamical generation of mass 
due to a magnetic field plays no role in the electroweak phase transition, because the field 
that had to be present in the early universe for the FDGM to be important was too large. 
Their claim was based in results obtained within QED4, where the critical field strength is 
~ 10 42 G. Here we have a critical field ~ 10 32 G. Although much smaller, it is still larger than 
10 24 G, which is the estimate generally predicted by primordial field generating mechanisms 
||13|| . We must point out, however, that the present result has been found within a toy model, 
so it just indicates the tendency of the theory when new interactions are switched on. We 
anticipate that in the electroweak model the critical field needed shall be much smaller than 
the one found here. 



The main conclusion of this paper is that Yukawa interactions enhance the dynamical 
generation of fermion bound states and masses in the presence of external magnetic fields. 




(43) 



For Yukawa coupling X y ~ 0.7 and gauge coupling g = |e we have 




(44) 



IV. CONCLUSIONS 
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It is natural to expect that other interactions will have similar consequences. Therefore, it 
is worth to extend our results to more realistic theories. In the context of the electroweak 
model the enhancement of the FDGM should be even more substantial due to the interaction 
richness of the theory. 

The generation of a fermion condensate in the presence of a constant magnetic field pro- 
vides an example of change in the symmetry properties of the vacuum due to external fields 
and may have very wide applications. The FDGM was originally discovered as a mechanism 
of catalysis of chiral symmetry breaking Q and as so, it has been rediscovered in many differ- 
ent theories, including condensed matter phenomena. We believe that in gauge theories with 
a richer set of interactions and symmetries the FDGM may be connected to the breaking 
of a gauge symmetry. In the context of the Standard Model our conjecture poses an inter- 
esting question: could primordial magnetic (hypermagnetic) fields induce gauge symmetry 
breaking in the early universe? If the answer is positive, it may led to several cosmological 
consequences. Work on this direction is in progress and will be publish elsewhere. 
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APPENDIX 



Let us find the solution of Z n that satisfies the equation 



z,i-P^-^l{\gH\) ' A ' 



(l + ZW-fe-o) 



(2tt) 4 q 2 (l + Z,,) 2 ^ - g|| ) 2 + m 2 ( 



(Al) 



Doing the variable change 
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k = (fa, fa) = (p 4 - q 4 ,p 3 - q 3 ) 



(A2) 



the integral in Eq. ( [All) can be written as 



d 4 q e"«x ( 1 + z ll h-(P ll -?„) 



(2tt) 4 g 2 (l + Z,,) 2 ^ - g|| ) 2 + m 2 ( 
<i 2 g_j_(i 2 A;|j 



(l + ^fr-*,, 



(27r)4 ^ + ^(Pii-M 2(1+ ^ )2( ^ )2+m2 ( fc n) 



(A3) 



We can change fa now to polar coordinates (k ,9^), and use Feynman integral formula 



dy- 



AB J -[yA+(l-y)B] 2 



to write the integral ( |7l3| ) as 

I = J wt~J d ^J dy J ddkH w e ~ z i 1 + z ") 
^ ' 

7 « H cos6» fc + 7 3 K n sin #fc 

X ■ 



(A4) 



where we have written the external parallel momentum in polar coordinates 



(A5) 



Pi = P4 = P|| COS^p 

Ps=P3 = P|| sin 6 P 
Using the following formulas, which are valid for a 2 7^ b 2 



dO- 



cos 9 



a sin 9 



+ 



a + bcos9) (a 2 — b 2 ) (a + b cos 9) a 2 — b 2 J a + bcos9 



dO- 



(A6) 



(A7) 



d9 



sin 9 



(a + bcos9f (a 2 -b 2 ) b (a 2 - b 2 ) 



In (a + b cos 9) 



(A8) 



d9 



2 arctg {a - b)t9 * 



a + b cos 9 



{a-b)tg\-^Jb 2 -a 2 
Vb 2 -a 2 111 {a-b)tgl+y/b 2 -a 2 ' 



In 



z'/ a 2 > b 2 
if a 2 < b 2 



(A9) 



and the definitions 
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~ 2 y 
a = yq± + 



2\gH\ 



(pf + + (1 - y) [(l + Z^iK^+m 2 



(A10) 



2p||K M 



(All) 



one can straightforwardly show that all the angle integrals are zero. Therefore, Eq. (Al) 
reduces to 



Z n7 -P„ = 



(A12) 



leading to the solution Z, =0. 



17 



REFERENCES 



[1] V. P. Gusynin, V. A. Miransky, and I. A. Shovkovy, Phys. Rev. Lett., 73, 3499 (1994); 
Phys.Lett. B 349, 477 (1995); Phys. Rev. D 52, 4718 (1995); ibid 4747; Nucl. Phys. B 
462,249 (1996). 

[2] C. N. Leung, Y. J. Ng, and A. W. Ackley, Phys. Rev. D 54, 4181 (1996). 
[3] D.-S Lee, C. N. Leung, and Y. J. Ng, Phys. Rev. D 55, 6504 (1997). 
[4] V. P. Gusynin and I. A. Shovkovy, Phys. Rev. D 56, 5251 (1997). 
[5] D. K. Hong, Y. Kim, and S-J. Sin, Phys. Rev. D 54, 7879 (1996). 
[6] D. K. Hong, Phys.Rev. D 57, 3759 (1998). 

[7] D. Ebert and V. Ch. Zhukovsky, Mod. Phys. Lett. A 12, 2567 (1997). 

[8] V. Elias, D. G. C. McKeon,V. A. Miransky, and I. A. Shovkovy, Phys. Rev D 54, 
7884 (1996); I. A. Shovkovy and V. M. Turkovsky, Phys. Lett. B 367, 213 (1996); D. 
M. Gitman, S. D. Odintsov, and Yu. I. Shil'nov, Phys. Rev D 54, 2968 (1996); A. V. 
Shpagin, |hep-ph/9611412j; I. A. Shushpanov and A. V. Smilga, Phys. Lett. B 402, 351 



(1997) ; E. Elizalde, Yu. I. Shil'nov, and V. V. Chitov, Class. Quant. Grav. 15, 735 

(1998) ; V. P. Gusynin, D. K. Hong, and I. A. Shovkovy, Phys. Rev D 57, 5230 (1998); 
V. A. Miransky, [hep-th/ 98051591 , V. P. Gusynin and A. V. Smilga, |hep-ph/9807486| . 



[9] E. J. Ferrer and V. de la Incera, |hep-th/9805T05| , to appear in Proc. of SILAFAE'98 



April 8-11, San Juan, Puerto Rico. 

[10] The FDGM may be important in the vicinity of magnetars. Magnetars are a rare type 
of star whose existence has been recently confirmed by astronomical observations of the 
largest batch of energy (in the form of gamma light) to arrive at Earth from a star 
beyond our solar system (see Kouveliotou, C, et al., International Astronomical Union 
Circular No. 7001 (Aug. 28) 1998). They are supposed to be responsible for the largest 



18 



magnetic field ever observed in our present universe, ~ 10 15 G. Such a field is about 100 
times larger than ordinary neutron star fields. 



11] K. Farakos and N. E. Mavromatos, cond-mat / 971028"B| ; K. Farakos, G. Koutsoumbas, 



and N. E. Mavromatos, Phys.Lett. B 431, 147 (1998); |cond-mat/9805402 



[12] G. W. Semenoff, I. A. Shovkovy, and L. C. R. Wijewardhana, |hep-ph/ 9803571 



[13] T. Vachaspati, Phys.Lett. B 265,258 (1991);K. Enquist and P. Olesen, ibid.319, 178 
(1993); B. Cheng and Olinto, Phys. Rev. D 50, 2421 (1994); G. Baym, D. Bddeker, 
and McLerran, ibid.53,662 (1996); P. Olesen. |hep-ph/9708320| , A. Brandenburg, K. En- 
quist, and P. Olesen, Phys. Rev. D 54, 1291 (1996); Phys.Lett. B 392, 395 (1997); K. 
Enquist JStro-ph/ 97073001 , G. Sigl, A. Olinto, and K Jedamzik, Phys. Rev. D 55, 4582 
(1997); M. Joyce and M. Shaposhnikov, Phys. Rev. Lett. 79,1193 (1997); M. Giovan- 
nini and M. Shaposhnikov,Phys. Rev. D 57, 2186 (1998); O. Tornkvist, Phys.Rev.D 
58:043501, (1998); |hep-ph/9801286t M. Hindmarsh and A. Everett, |astro-ph/ 9 708001 ; 



D. Grasso and A. Riotto, Phys. Lett. B 418, 258 (1998); J. Ahonen and K. Enquist, 
Phys. Rev. D 57. 664 (1998). 



[14] J. Ambjorn and P. Olesen, [hep-ph/9304220 



[15] E. J. Ferrer and V. de la Incera, Phys. Rev D 58, 065008 (1998). 
[16] J. Schwinger, Phys. Rev. 82, 664 (1951) 

[17] V. I. Ritus in Issues in Intense-Field Quantum Electrodynamics, ed. V. L. Ginzburg 
(Nova Science, Commack, 1987). 

[18] Handbook of Mathematical Functions, eds. M. Abramowitz and I. A. Stegun (Dover, 
New York, 196). 

[19] A straighforward calculation can convince us that the fermion condensate in the model 
considered in this work has the same expression as the fermion condensate found in 



19 



ref. [3], that is, (ipip^ = iTr ^G(x, x)\ ~ —^^-m In (-^Pj • The difference now is that 
the present mass is given by Eq. (fn|), so the condensate is enhanced by the Yukawa 
coupling. 



20 



